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Abstract 

The propagation of optical pulses in two types of fibers with randomly varying 
dispersion is investigated. The first type refers to a uniform fiber dispersion su- 
perimposed by random modulations with a zero mean. The second type is the 
dispersion-managed fiber line with fluctuating parameters of the dispersion map. 
Application of the mean field method leads to the nonlinear Schrodinger equation 
(NLSE) with a dissipation term, expressed by a 4th order derivative of the wave 
envelope. The prediction of the mean field approach regarding the decay rate of a 
soliton is compared with that of the perturbation theory based on the Inverse Scat- 
tering Transform (1ST). A good agreement between these two approaches is found. 
Possible ways of compensation of the radiative decay of solitons using the linear and 
nonlinear amplification are explored. Corresponding mean field equation coincides 
with the complex Swift-Hohenberg equation. The condition for the autosolitonic 
regime in propagation of optical pulses along a fiber line with fluctuating disper- 
sion is derived and the existence of autosoliton (dissipative soliton ) is confirmed 
by direct numerical simulation of the stochastic NLSE. The dynamics of solitons 
in optical communication systems with random dispersion-management is further 
studied applying the variational principle to the mean field NLSE, which results in 
a system of ODE's for soliton parameters. Extensive numerical simulations of the 
stochastic NLSE, mean field equation and corresponding set of ODE's are performed 
to verify the predictions of the developed theory. 
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1 Introduction 



Propagation of optical pulses in fibers with non-uniform dispersion remains 
to be a field of intensive research. The relevant studies are motivated by their 
application to modern optical communication systems. Random fiuctuations 
of the dispersion coefficient is inherent to existing optical fibers, which dete- 
riorates the performance of communication lines based on both standard and 
dispersion-managed solitons [1]. Therefore, the proper compensation of the 
pulse distortion caused by random fiuctuations of the fiber dispersion is of 
vital importance, especially for long haul optical communication systems. 

The previous investigations were concerned with the analysis of the modula- 
tional instability (MI) of electromagnetic waves in fibers with random disper- 
sion. It was shown that the randomness can reduce the MI gain for anomalous 
dispersion region and generate the instability for the whole spectrum of mod- 
ulations in the normal dispersion region [2,3,4]. Recently the performance of 
dispersion-managed soliton system with random variations of span length and 
span path-average dispersions has been studied numerically in Ref. [5]. The 
comprehensive analysis of the dynamics of solitons implies the solution of the 
nonlinear Schrodinger equation (NLSE) with randomly varying coefficients, 
which is a rather complex problem. However, in the particular case, when a 
uniform fiber dispersion is superimposed by weak random modulations, the 
perturbation theory based on the Inverse Scattering Transform (1ST) can be 
applied. In the framework of this approach the decay law for optical solitons 
was derived which describes their radiative damping due to emission of linear 
waves [6,9]. The case of dispersion-managed (DM) solitons is more compli- 
cated because here one has the superposition of the random fluctuations and 
strong periodic modulations of the dispersion. In recent works [10,11,12,13] the 
variational approach was employed to analyze the dynamics of optical pulses 
in systems with random dispersion- management. It was shown that the DM 
soliton is subject to disintegration due to random variations of the dispersion 
map parameters. The averaged equation in the frequency domain was derived 
in [14,15,16], where an interesting idea of pinning (compensation of the ac- 
cumulated effect of fiuctuations) is proposed. Numerical simulations of the 
dynamics of DM solitons in systems with random dispersion-management are 
presented in [10,11,19,20]. Analytical description of the DM soliton dynamics 
in fiuctuating media remains to be an open problem, since the radiative effects 
cannot be taken into account in the frame of the variational approach. How- 
ever, there exists a limit pertaining to strong dispersion-management, where 
the NLSE with periodic coefficients is nearly integrable [21]. The development 
of the stochastic perturbation theory similar to the 1ST approach may be 
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successful in this limit. This problem deserves a separate consideration. 

In this study we apply the mean field method (MFM) for the analysis of 
propagation of optical pulses in fibers with random dispersion, including the 
DM case. The advantage of this method consists in that, it doesn't depend on 
the fact whether the basic deterministic equation is integrable or not. 

The paper is organized as follows. In section II we derive the mean field equa- 
tion for the case of uniform dispersion perturbed by random modulations and 
compare the predictions of this approach with the corresponding result of the 
1ST based perturbation theory. The condition for the existence of autosoliton 
in random media is derived. The section III is devoted to analysis of the DM 
soliton dynamics in the framework of MFM. We show that the variational ap- 
proach applied to the mean field equation with a frequency dependent damping 
leads to a system of ODE's for DM soliton parameters. And in the last section 
IV, we briefly summarize the main results of this study. 



2 Mean field method and 1ST approach 

The propagation of optical pulses in a fiber with uniform dispersion superim- 
posed by random modulations was considered in the framework of the pertur- 
bation theory based on the 1ST [6]. Assuming the perturbation to be of the 
form R = e{z)uu, where e{z) is a weak random process, the radiative decay of 
a soliton was calculated. The decay law for the pulse amplitude was derived 
from the energy conservation and the dynamic balance between the radiative 
component and localized mode. It was shown that the amplitude of the soliton 
decays with the propagation distance as for z < l/o"^, where o"^ is the 

noise level. In addition, the decay rate was revealed to be highly dependent 
on the pulse duration ~ t^^, i.e. the infiuence of the randomness is superior 
on short pulses, particularly in the femtosecond range. 

The extension of this approach for more general cases (e.g. dispersion-managed 
solitons) is in general not possible. Therefore, it is of interest to apply another 
approach which is independent on the fact whether the underlying determin- 
istic equation is integrable or not. One of such possibilities is provided by 
the mean field method [22]. Despite the well known troubles of this approach 
for random nonlinear wave equations [23,24], one can obtain reasonably well 
description in the case of weak fiuctuations, and/or small distances of propa- 
gation. 

In this section we apply the MFM and the perturbation theory based on 
the 1ST to the problem of optical pulse propagation in fibers with random 
dispersion. By comparing the predictions of these two approaches we reveal 
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the limits of validity of the MFM. 



Optical pulse propagation in fibers is well described by the NLSE 

iu^-l — ^-^Uft + \u\'^u — i5u + iijL\u\'^u. (1) 



Here 6, jj, are the coefficients of linear and nonlinear amplification (damping) 
respectively. In our model these terms are supposed to compensate the damp- 
ing originated from the randomness of the fiber dispersion. Below we consider 
the case, when the dispersion coefficient d{z) is the sum of the constant and 
random parts, i.e. d{z) = 1 + e{z) with 

< e 0, < e{z)e{z') >= B{z - z', iX-^o ^ <j''S{z - z'). (2) 



The noise is assumed to be small compared to the constant part of the dis- 
persion. 

To derive the equation averaged over fluctuations we apply the mean fleld 
method, representing the field as consisting of the mean value and small fiuc- 
tuating part u —< u > +5u,< Su >— 0,5u <S< u >, where < it > is a 
slowly varying mean field. According to this method we can use the decou- 
pling < \u\'^u >^ I < M > p-u. This corresponds to neglecting by fluctuations 
of the nonlinearity, and the approximation is valid for propagation distances 
z « In typical optical experiments this requirement is satisfled. In gen- 

eral the decoupling procedure is inaccurate since the scattered field 5u grows 
with propagation distance and becomes comparable with < u >, that violates 
the assumption 6u <^< u >. In order to decompose the mean < e{z)utt > we 
use the Furutsu-Novikov formula 

< e{z)F{u) >= [b{z - z') < ^^TT > d^'- (3) 



Of particular interest is the case of white noise fiuctuations, when B{x) ~ 
a'^S{x). Also the causality principle will be useful 

^ = 0, a z'< 0, z' > z. (4) 



Then integrating Eq.(l) over z' from to z and taking variational derivative 
(with the causality principle in mind) we obtain 

< e{z)uu{z, t)>^ — <u >tttt ■ (5) 
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In the result we get the equation for averaged pulse profile with a fourth order 
dissipative term (for simplicity we dropped averaging symbol in < u{z, t) >) 



iUz H — Utt + \u\'^u = —i^Utttt + 'i^i^ + ilJ'\u\'^u, (6) 
2 



where 7 = a being the noise strength. 

As it is apparent from this equation, the effect of random dispersion on the 
pulse evolution is described as the pulse propagation in a uniform medium 
with the frequency depending dissipation. Formally this equation coincides with 
the complex Swift- Hohenberg equation. Thus we can expect the existence of 
dissipative solitons in nonlinear media with fiuctuating dispersion. 

At first we examine the decay rate of a pulse under the action of the 4th order 
dispersion coefficient in Eq.(6). For the weak noise case this term is small 
and we can apply the perturbation theory. The single soliton propagation is 
conveniently described using the equations for the energy N — J \u\^dt and 



momentum P — Im(/ uludt). 



dN 
dz 

dP_ 
dz 



:2 j [5\u\'' + ix\u\''-^\uu\'']dt, (7) 
—00 

00 

2Im j [5u + iJi\u\^u + ^Uttttl'^t^t. (8) 



We look for the single soliton solution in the form 

u{z, t) = r/sech[r;(t + l]^)]e-^"*+^(^'-"')^/'. (9) 



Substituting Eq.(9) into Eqs. (7), (8) we obtain the system of equations for the 
soliton amplitude and velocity 



^=2„|* - ^(lg,f + 2,fa^ + Si') + Hm'l, (10) 

Parameters 5, /i are absent in the second equation, since the soliton velocity is 
not affected by such type of perturbations. For ^2 = one can find the decay 
law of the soliton due to the 4th order dispersive dissipation, which follows 
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from Eq.(lO) 



Now we can compare the prediction of the MFM with that of the perturbation 
theory based on the 1ST [6] and thus find the hmits of apphcabihty of the MFM 
(within the hmits of the error of the 1ST method). In the frame of 1ST based 
perturbation approach, the decrease of the sohton amphtude due to radiation 
of random waves under fluctuations of dispersion can be estimated from the 
balance equation for the norm N — J \u\'^dt, which is the integral of motion. 
The norm can be represented as 

oo 

N^2r) + - f Ln\a{X)\-'^d\, (13) 

TT J 

— oo 



where A = 2A; is the spectral parameter, k is the wavenumber, a{X),b{\) 
are the Jost coefficients of the Zakharov-Shabat linear spectral problem and 
|ap + |6p = 1. For the weak noise case Ln(|a|~^) ^ |&(A)p. The balance 
equation can be obtained by differentiation of the norm with respect to z 



d<N> 
dz 



oo oo 

-i/Q(A)dA-i/(|6a-^dA + + (14) 



where 



rlh*(\) 

g(A) = 2Re< 6(A)— ^ >, 

dz 



is the spectral power of emitted radiation by the soliton. The Jost coefficient 
6(A) can be calculated using the perturbation theory based on the 1ST [7,8]. 
It should be noted that the total wave is given by 

u{z,t) = Us{z,t) +Uc{z,t), (15) 



where 



in J a (A + irjY 



z/2 cxp(2/(2((f-f,) + oJ) 



ITT 



cosh (2;) 2 



7 b* 1 

J a*^ 'X-ir] 



(16) 
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Here Ug is the soliton part of the total wave. The first component of u,. rep- 
resents the emitted wavepacket. The second component is due to interaction 
between the sohton and the emitted wavepacket and it is appreciable only 
nearby the sohton. We will be interested in the decay of a soliton under emis- 
sion far from the soliton, assuming that this emission is lost and has no back 
action on the soliton. 



The equation for the amplitude is 

-;- = 7, e = — cr . (17) 

dz l + 5eV 15 ^ ' 



The solution reduces to the algebraic equation 

T] ezjf = r]Q. (18) 



Since r)— < r) >~ e we can replace < rj^ >!^< rj >^. 

For distances of propagation z < l/cr^ one can neglect the second term in the 
denominator of Eq.(17) and find the decay law for the amplitude 

< "^^'^ (1 + 1^X^.)V4 - ^^0(1 - Tb^'^'^oZ + 0{a^)). (19) 



This equation gives almost the same decay rate as Eq.(12) (see [18]). 

Thus, the 1ST approach conforms with the picture of soliton propagation in 
uniform media with the effective frequency dependent (~ uj'^) losses. Also the 
1ST and MFM predict the strong ~ l/t^ dependence of the decay rate on 
the pulse duration. The comparison shows that the mean field theory un- 
derestimates the soliton decay rate relative to the 1ST as given by the ratio 
P = r MFM 1ST = 0.875. This discrepancy between numerical constants in 
Eq.(12) and Eq.(19) is due to the approximate character of the mean field 
equation, where we have neglected by the renormalization of the nonlinear 
term. 

In Fig.l we report the comparison of predictions of different models regarding 
the decay law of the pulse (decreasing of its amplitude) in the course of prop- 
agation along a fiber line with randomly varying dispersion. It is seen that all 
the models well describe the initial stage of the radiative damping. The MP 
NLS Eq.(6) and decay law Eq.(12) better describe the damping at distances 
z > l/cr^, compared to the algebraic Eq.(18). The curve labeled stochastic 
NLS is obtained by the averaging over 70 reahzations. 
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Fig. 1. Decay law for the pulse amplitude according to different models. The initial 
pulse has the form u{0,t) = riQsech.{riQt) with rjQ = 1. The fiber dispersion randomly 
varies around do = 1 with a mean square fluctuations a = 0.3. 

Numerical integration of the stochastic NLS Eq.(l) and MF NLS Eq.(6) are 
performed by the split-step fast Fourier transform. Absorption on the domain 
boundaries is employed to imitate the infinite length. The system ODE's 
Eq.(lO), Eq.(ll) and Eq.(28) are solved using the procedure D0PRI8 [25], 
which is based on the Runge-Kutta scheme with adaptive stepsize control. 

Inspecting the equation for the field momentum P one can see that Qz = 0. 
This is due to the fact that the numbers of quanta emitted by the soliton in 
forward and backward directions are equal. Therefore the predictions of the 
MFM and 1ST for Q coincide only at the fixed point, where = 0. The ratio 
of amplitudes at the fixed point is r = timfm/vist ~ 0.97 for the case of linear 
amplification, and r = 0.93 for nonlinear amplification. 

The decay of a soliton can be compensated by linear and/or nonlinear amplifi- 
cations. This results in formation of an autosoliton. In our case, when the pulse 
power loss due to 4th order dispersion effect i'-fUutt is compensated by the lin- 
ear i6u and/or nonlinear ifi\u\'^u amplification in Eq.(l), such an autosoliton 
can be created. A distinctive property of these solitons is that, they recover 
the stable waveform when deformed. Fig.2 demonstrates such an aspect of 
autosolitons, when the power dissipation due to 4th order dispersion term is 
compensated by linear /nonlinear amplification. Note that the linear amplifi- 
cation gives rise to more pronounced oscillations around the fixed point during 
transient period compared to nonlinear amplification. At longer propagation 
distances the oscillations are dumped out and a stable dissipative soliton is 
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formed. 



1.6 - 

^5 _ (a) a = 0.3, 8 = 0.01, n = 0. 




25 50 75 100 125 150 175 200 
Distance (z) 



1.6 - 

15 . (b) a = 0.3, 8 = 0, n = 0.015. 
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Distance (z) 



Fig. 2. Autosoliton acquires its stable form at some propagation distance when its 
amplitude initially assigned the value greater or lower than corresponding fixed 
point value r] = 1.0. Obtained by numerical solution of the MF NLS Eq.(6) for 
a = 0.3, 7 = = 0.0225, with: (a) Hnear amplification S = 0.01, = 0.0, and 
(b) nonlinear amplification 6 = 0.0, fi = 0.015. 



The amplitude of the autosoliton at the fixed point can be found from the 
balance equation (14) of the 1ST 



Two types of fixed points exist. The first type rjist is defined by the competition 
between the dissipation (induced by the randomness of the fiber dispersion) 
and the hnear/nonhnear amplifications. 

For 5 — 0, II > (the nonhnear amphfication dominates) we find 

vL = % (21) 



and similarly for = 0, 5 > (linear amplification dominates) 



2 



The second type of fixed points is defined by the competition between the 
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linear dissipation and nonlinear amplification. 



?72st 



3\6\ 9|(5|V ^ 
2n 20/^3 



(23) 



i.e. the fluctuations reduce the value of amplitude in this case. 

The fixed points for the soliton amplitude (at a — 0.3, S — 0.01, = 0.015) as 

predicted by Eqs.(21), (22) are as follows : rjut = (S/i/a^)^/^ ~ 0.91, and rjist — 
(155/2ct2)1/4 ~ 0.96. The fixed points of ODE (10) for the same parameters 
give T] = {AOfi/la^/^ ~ 0.976, r] = {606/7aY/^ ~ 0.988, which is the better 
approximation to PDE simulation of Fig.2, when we assume that the effect of 
noise a is accounted by the 4th order dispersive dissipation 7 = (T^/4. 

Simulations of the stochastic NLS Eq.(l), is presented in Fig. 3. Compensa- 
tion of the damping (due to the randomness of the dispersion) by a linear 
amplification gives rise to a stable soliton. To verify the property of autolitons 
observed in Fig.2, we assigned the amplitude greater and lower values, com- 
pared to the stationary one Ugt — 1.06. As expected, the soliton adjusts itself 
to the stationary amplitude at some propagation distance. Similar behavior is 
observed when the combination of linear amplification and nonlinear damping 
is applied (Fig. 3b). However, it should be pointed out that the stable pulse 
propagation is limited by the growth of the zero mode under amplification. 
When the initial pulse amplitude is close to this solution the instability starts 
to at distances z r^l/ji (see f.e.[26]). 




Distance (z) Distance (z) 



Fig. 3. Autosoliton from stochastic NLSE simulations, (a) Evolution of the soliton's 
amplitude towards the stationary value as obtained by numerical solution of the 
stochastic Eq.(l) (solid lines), and as predicted by the ODE model Eq.(lO) (dashed 
lines). Dash-dot line is the stationary value of the amplitude, (b) Autosoliton re- 
sulted from the combined effect of a linear amplification and nonlinear dissipation. 
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The phase trajectories according to Eq.(lO) and Eq.(ll) is shown in Fig.4. 
The fixed point appears to be of sink type. 




Fig. 4. The phase trajectories corresponding to ODE model Eq.(lO) and Eq.(ll). 
a = 0.1, (5 = 0.0015, n = 0. 

Thus, summarizing the results of this section we note, that the optical pulse 
propagation in fibers with random dispersion can be described as propagation 
in a uniform fiber with the effective frequency dependent damping. The effect 
of randomness of the fiber dispersion is accounted by a term with 4th order 
derivative of the wave envelope. Application of this idea to dispersion-managed 
solitons will be considered in the next section. 



3 Dispersion-managed soliton in a fiber with random dispersion 

The dispersion-managed optical communication line consists of fiber sections 
with alternating anomalous and normal dispersion coefficients [27]. Now the 
function d{z) is the sum of periodicaUy varying dispersion dQ{z) and the ran- 
dom part di{z) i.e d{z) = do{z) + di{z). The usual scale of the DM map is ~ 
100 km, so we can consider the random modulations of the dispersion in the 
white noise limit with 5{z) correlation. Due to strong variations of the disper- 
sion within the unit cell (the strong dispersion-management regime), locally 
one has the linear dynamics. The effect of nonlinearity is small in the scale of 
a unit cell, while accumulated on longer distances, can lead to formation of a 
breathing DM soliton. Therefore we can expect that the inffuence of the ran- 
dom modulations of dispersion will result in the frequency dependent damping 
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in the form calculated in the previous section and the nonlinear corrections to 
this damping will be small. 

Performing the averaging over fluctuations yields the following equation 

iuz + ^^^^utt + \u\'^u = -i^utttt + iSu + iiJ,\u\^u. (24) 



In the case of strong dispersion-management (when d{z) — {l/e)d{z/e), e <S 1, 
the averaging procedure can be justified [28]. Introducing the variable ( — z/e 
we obtain the equation 

iu^ + — T^-Mtt + — i6 — iiJ,\u\'^)u = 0, (25) 



where d{z) = (io(C) + di{C)i di ~ e^^^. Let us consider the case when di ~ ^/e 
and thus dl ~ e. One can decompose the field as u =< u > +5u, 6u ~ y^, 
and obtain the following estimates for the correction l^wp, 6u'^ ~ e. Thus, the 
nonlinear corrections from fluctuations to the averaged equation are of order 

and can be neglected in comparison with terms < di{Quu >~ e\ < u > 

<u >~ e. 

The variational approach is proved to be effective for exploring the dispersion- 
managed soliton [29,30]. In the presence of nonconservativc terms in Eq.(24) 
we can use the modified variational equations. The equations for the pulse 
parameters are [31,32] 

dL d dL _ 7c^^"*_^^^^ C2g) 
drji dz drji J drji 

— oo 



Here L is the averaged Lagrangian L — J L{z,t)dt. Employing the Gaussian 
anzats for the waveform 

uiz, t) = A{z) exp[-^-^ + + tcPiz)]. (27) 

we obtain the following system of ODE's from the Eq.(26) 



{A% = -d,{z)A% - \lA\^, + 26^ - h'a'\ + 2M^ + 

= ^ - d,{z)}? - - 127K:^ + ^')- (28) 



or \/Za^ a" 
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When ^ — 5 — IjL — 0, we reproduce the variational equations for a DM sohton 
[29]. 



In order to verify the vahdity of the mean field equation (24) and correspond- 
ing variational equations (28) for description of the DM sohton dynamics 
we compared the result of their numerical solution with that of the original 
stochastic equation (1), averaged over 50 realisations of random paths. As an 
initial condition we take the waveform Eq.(27). Parameters of the dispersion 
map are: Di — 30, D2 = —28, zi = Z2 = 0.078, Zd = zi + Z2. Fluctua- 
tions of the fiber dispersion coefficient is supposed to be of type Eq.(2), and 
d{z) = do{z)[l + e{z)]. Satisfactory agreement of these data is displayed in 
Fig. 5. A characteristic feature of the pulse decay process due to the 4th order 



a = 0.05 




Distance (z/z) 

d 



Fig. 5. (a) Comparison of pulse decay rates obtained by numerical integration of the 
stochastic NLS Eq.(l) (broken line), the averaged equation Eq.(24) (solid line), and 
system of ODE's (28) (dashed line). DM map parameters and initial conditions arc 
Di = 30, D2 = -28, zi = Z2 = 0.078, = 0.156, A{0) = 1.52, a(0) = 1.059, 6(0) = 0, 
and7 = c7V4, a = 0.05 ^0.1, 6 = i^ = 0. 



dispersion term is that, the pulse amplitude rapidly decreases at initial stage, 
after that almost linear law is followed. This behavior can be understood from 
the fact, that Uutt term depends on the pulse shape: its value is bigger when 
the pulse is more sharp in the initial stage, and becomes small when the pulse 
is dumped and broadened. 

Now wc try to compensate the pulse damping due to randomness of the fiber 
dispersion (which is accounted in the mean field equation as a 4th order deriva- 
tive term), by linear and/or nonlinear amplification. The Fig. 6 demonstrates 
that this is indeed achievable. 
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= 0.1, 6 = 0.06, n = 0. (a) 
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Fig. 6. (a) Compensation of the pulse decay by means of linear amplification gives 
rise to a stable DM soliton: broken line - stochastic NLS Eq.(l), dashed line - ODE 
(28), solid line - mean field Eq.(24). (b) Combined action of linear and nonlinear 
amplifications on a longer distance DM soliton propagation. The inset shows the 
dispersion map (unperturbed -solid line, purturbed - dashed line). DM map param- 
eters and initial conditions are similar to those of figure 5. 

An important kind of randomness inherent to real dispersion-managed systems 
is related to the fluctuations of the lengths of fiber spans Zi, Z2. In Fig. 7 we 
illustrate that the decay of a DM soliton due to the randomness of fiber span 
lengths can be compensated by a suitable linear amplification. Comparison 
with the non-amphfied case {6 — 0) shows the effective stabihzation of the 
DM soliton. 




10 20 
Distance (z/z .) 



10 20 30 
Distance (z/z ) 



Fig. 7. (a) Decay of the amplitude of a DM soliton at different strengths {a = 0^0.2) 
of fluctuations of the fiber span lengths zi ■ (1 + e(^)), Z2 ■ {1 + £{z)) according to 
Eq.(l) with 6 = f^ = 0. (h) Compensation of the pulse decay (at a = 0.2) by means 
of a linear amplification {S = 0.075, = 0) gives rise to a stable DM soliton. The 
inset shows the dispersion map (unperturbed -solid line, purturbed - dashed line). 



Now let us analyze the fixed points of the system. These points correspond to 
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the stable propagation of a DM soliton under fluctuations of the dispersion 
map parameters. Equating the left hand sides of Eqs.(28) to zero and assuming 
6 ~ 5, /X, 7 -C 1, we obtain 



(29) 




(30) 



1 (biiAl 67 



(31) 



'° 2(io ^ 2 V2 4 



According to Eq.(29), the stable DM soliton propagation in optical flbers 
with random dispersion is possible only in the region of anomalous dispersion 



4 Conclusion 

We have studied analytically and numerically the propagation of optical soli- 
tons in fibers with random dispersion. The cases when the unperturbed disper- 
sion is uniform and periodically modulated (DM) are considered. It is shown 
that the optical pulse dynamics can be successfully described in the framework 
of the mean field method. The limits of validity of the MFM is revealed by 
comparison with the results of the perturbation theory, based on the 1ST. It 
is shown that when the linear and nonlinear gain/damping are included, cor- 
responding mean field equation coincides with the Swift-Hohenberg equation. 
Analysis revealed the existence of dissipative solitons in fibers with fiuctuating 
dispersion. The system of variational equations for parameters of the DM soli- 
ton has been derived, which takes into account the nonconservative effects due 
to the randomness of the fiber dispersion. We have analyzed the fixed points 
of this system and found conditions for the stable propagation of DM soli- 
tons in a fiber with gain and fiuctuating dispersion. The analytical results are 
confirmed by direct numerical simulations of the full stochastic NLS equation. 
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